Initial wave packets and the various power-law decreases of 
scattered wave packets at long times 
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The long time behavior of scattered wave packets tp(x,t) from a finite-range potential is investi- 
gated, by assuming ip(x,t) to be initially located outside the potential, ft is then shown that ip(x,t) 
can asymptotically decrease in the various power laws at long time, according to its initial charac- 
teristics at small momentum. As an application, we consider the square-barrier potential system 
and demonstrate that ip(x,t) exhibits the asymptotic behavior t~ 3 ^ 2 , while another behavior like 
t~ 5 / 2 can also appear for another ip(x,t). 

PACS numbers: 03.65.-w, 03.65.Nk 



A long time deviation from the exponential decay law 
is predicted in unstable quantum systems such as an a- 
decaying nucleus 0. These systems are often modeled 
in the systems of a particle in a short range potential. In 
this approach, a nonexponential decay is generally found 
to follow a particular power law, £~ 3 / 2 (see, e.g., Refs. 
00,0, On the contrary to the theoretical results, 
an experimental evidence still has not been discovered [6j , 
which requires us to reexamine such a nonexponentially 
decaying behavior from every aspect. As far as the au- 
thor knows, conventional studies are developed without 
taking into account of the characteristics of the initial 
wave packets. A consideration for them may lead to a 
discovery of other aspects of the subjects. Indeed, as 
for the free particle system, the various power decreases 
of wave packets, not restricted to i -3 / 2 , have recently 
been demonstrated and shown to be characterized by 
the behavior of initial wave packets at small momentum 
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In this paper, we analytically prove that the wave 
packet ip(x,t) scattered from a finite-range potential can 
asymptotically behave like t~^ 2 at long times, where 
j = 1,2,.... By making an assumption that the ip{x,t) 
is initially located outside the potential, these various 
power behaviors can be characterized by the initial char- 
acteristics of the ijj(x,t) at small momentum, like for the 
free particle case [13 • The validity of our analysis 
is numerically confirmed, through the application to the 
system with a square barrier potential. 



For a one-dimensional system with a potential V{x), 
the Hamiltonian H is defined as H = Hq + V, where 
H = —{h 2 /2M)d 2 /dx 2 being the free Hamiltonian. For 
simplicity, we use the unit such that h — 1 and 2M = 1 
throughout the paper. Although we confine ourselves to 
the one dimensional case, the following discussion can be 
extended to that in an arbitrary dimension in principle. 
The potential V(x) is assumed to have a finite-range, i.e., 



V(x) = for |z| > R, 



(1) 



for a positive number R. We only consider the systems 
without bound states. This restriction however is easily 
moderated. By using the eigenfunction expansion, the 
wave packet ip(x, t), evolving from the initial state (wave 
packet) ip{ x )> is expressed as 



1>(x,t) = (e-« u il>){x) = 



<p(x,k)i/>(k)dk, (2) 



where the i>{k), determining the initial energy- 
distribution, is defined by 



tp(k) = / (p*(y,k)ip(y)dy. 



(3) 



The functions (f(x, k) (k € M\{0}) are stationary scatter- 
ing solutions of the time-independent Schrodinger equa- 
tion, 



[H + V]<p(x,k) = k 2 ip(x,k). 



(4) 



More precisely, they must satisfy the Lippmann- 
Schwinger equations in one dimension \H\ , 



tp(x,k) 



1 



e *m\*-y\ V {y) tpfak) dy. 



/2ft 2i|*| 

(5) 

Either solution of the above equation with (+) or (— ) sign 
is accepted for ip(x,k) in Eq. (J2J. Notice that, being in- 
tegral equations, the Lippmann-Schwinger equations al- 
ready incorporate the boundary conditions, unlike Eq. 
(QJ. We here choose the equation with (+) sign. Then, 
tp(x, k) for positive k and for negative k are solutions of 
Eq. |@}, to the cases of an incident plane wave from the 
left and from the right of the potential, respectively. 

We first derive an asymptotic expansion of ip(x,t) at 
long times. The integral in Eq. @ is changed to a 
Fourier- integral form over the energy variable E = k 2 , 
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(6) 



2 



where T±(x,E) = ip(x,±E 1 / 2 )ip(±E 1 / 2 ). We may de- 
compose T±{x,E) into the following forms, 



as E — > 0, respectively. Then, substituting them into Eq. 
(|1U[) leads to an expected result that 



F±(x,E)=±E 1 / 2 ± (x,E)+£ ± (x,E), (7) ^ 



where 



°° p{r+s-l)/2 _ 

o±(x,E)= — n — «*,±o)v> (s) (±o), 



r+s— odd 



£±(x,£) 



E • 

r+s— even 



7-1 C I 



(8) 



9^(:r,±0)^ (s) (±0), (9) 



as -B — * 0. Both y>(x, fe) and tp{k) are assumed to be 
differentiable with respect to k without the origin. Fur- 
thermore, we have used the notations that dif>{x, ±0) = 
limfc^±o d n (p(x, k)/dk n , and for any function of k, say, 
f(k), f^(±0) = \im k ^ ±0 d n f(k)/dk n . Note that by set- 
ting Eq. (JTJ into ©, 0±(x,E) makes no singularity at 
E = while E~ 1 / 2 £±(x,E) does not necessarily. Then, 
in taking such a singularity into account, the asymptotic 
form of the Fourier integral © may read formally [l2| 



ip(x,t) 



1 1 

TTTTrr a&LOJx,0) 

2 ^ UtV + 1 ^ E a\ , J 

j=0 y ' 



<7 = ± 



i - r(j + i/2) ^ , 



(10) 



as t 



oo, where liniE-xx, d E 0±(x, E) = and 
lim^^oo ff E £±{x, E) = were assumed. It is noted 
that we could also obtain the expansion 1)1 OJI by different 
methods (see, Refs. 00 [10 and references therein). 

As same as the results in Refs. 011 IE EH, we 
can expect from Eq. IjlUI) that by specifying the low- 
energy behaviors of both 0±(x,E) and £±(x,E) ap- 
propriately, tp(x,t) can asymptotically show the various 
power-decreases. This situation may be realized by such 
an initial wave packet ^>{x) that satisfies 



^(™)(±o) = for n = 0,l,...,m-l, 



(11) 



with a certain integer m. In fact, let us first consider the 
case of m being an even number given by m = 2m. In 
this case, 0±(x,E) and £±(x,E) read 



0±(x,E) 



e±(x,E) 



fim/2 



9 fe (^(x,±0)V (m) (±0) + 0(£ (m+2)/2 ), 

(12) 



fim/2 



tp(x, ±0)i/; (m) (±0) + <3(£ (m+2)/2 ), 

(13) 



cr=± 



^kr^(x,0) 



r(m+l/2) 

m!( . i)w+1/2 c E M^uj 



+ o(t- m - 3/2 ), 



(14) 

where d E l O±(x,0) = m\d k ip{x, ±0)^ ( " i) (±0)/ to! and 
5f£±(x,0) = m!(p(x,±0)V' (m) (±0)/TO!- On the other- 
hand, if m is an odd integer m = 2m — 1, one see that 

p(m-l)/2 

0±(x,E) = : (^(x,±0)V (m) (±0) + O(^ (m+1)/2 ), 

m! 

(15) 



E (m+l)/2 

£fc(z,E) = (m + l)^(x,±0)V (m) (±0) 

(m + 1)! L 

+<^(x,±0)^ (m+1) (±0)l + 0(£ (m+3)/2 ), 

(16) 

as E — ► 0. Inserting them into Eq. I|1(J|) again, we obtain 



V>(x, i) 



1 



CT=± 

r(m+l/2) , ' 

- !(it p +1/2 ^M^0) 



0(f 



(17) 

where d^~ 1 ± {x^) = (m - l)\<p(x, ±0)^ (±0) /ml 
and df£±(x,0) = m![(m + l)fljk¥>(a, ±0)^ m )(±0) + 
^(x,±0)^ (m+1) (±0)]/(m + 1)!. 

We have to give such initial wave packets ^>(x) satisfy- 
ing the formal condition (|ll|l . Then, it seems practically 
advantageous to rewrite this condition in terms of the 
initial wave packet ip(k) in momentum (or Hq) represen- 
tation. This can be achieved by assuming the ip(x) to be 
located in the left of the scattering potential V(x), i.e., 



tp(x) = for x > —R. 



(18) 



This assumption will however be relaxed to that ip( x ) = 
for \x\ < R, in the discussions below. From the assump- 
tion Q18)). the i>(k) is expressed by the integral over the 
truncated interval (-co, — R), 



m = 



-iky 



ip(y)dy = 



-R e -iky 



/27T 



1>(y)dy. (19) 



Meanwhile, the assumption implies that <p(x, k) for 
x < — R is written by a superposition of plane waves 



<p(x, k) = [g+{k)e^ x + g.(k)e- l ^ x }/V2^. 



(20) 



Since we adopt Eq. J5J with (+) sign, we see that g+(k) = 
1 for k > and for k < 0. Substituting Eqs. JTjjJ and 
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(12011 into ®, we can obtain a desirable expression for 
as a linear combination of -0^(±O)'s, 

^(<r0) = {<rl) n 5+,„ ^ (n) (+0) 

+E(")(-- 1 )'3* (n ~Vo) ^(-o), ( 21 ) 

where aO stands for the limit symbol +0 or —0 for a = + 
or — , respectively, and <5 +i(T denotes Kronecker's delta. 
From this, one might expect that the condition at 
low energy, which determines the asymptotic forms i|14|) 
and 1|17|) . is characterized by the condition at small mo- 
mentum [T^j . 

$ {l) (±0)=0 for / = 0,1,..., m- 1. (22) 

This condition indeed implies the condition How- 
ever the converse does not necessarily hold. This incom- 
patibility comes from the actual behavior of g^~ l \±0ys. 

Let us now demonstrate the asymptotic formulas 114|) 
and 1)17(1 with Eq. (|21|l , by applying it to the system with 
the square barrier potential V(x) given by Eq. Q and 



V(x) = V a for |x| < R, 



(23) 



where Vq (> 0) is the height of the potential barrier. For 
this system, g-(k) in Eq. (|2"0")) is given by (see, e.g., 0) 



ff-(*0 



2ikp 
9{-k) 



-g(k) sinh 2pR for < k < k b 



for - k b < k < 



where k b = Vq 1/2 , p = [k b 2 - k 2 ] 1 ' 2 , and 



g{k) 



cosh2pi? 



2ikp 



sinh 2pi? 



-2ikR 



(24) 



(25) 



The g(fc) has the remarkable properties: \inik_ t ± g(k) = 
and limfc^ (T o.9(' 7 fc)/fc — o2j(ik b sinh 2k b R), to give 



0_(+O) = -l, ff _(-0) = 0. 



(26) 



This implies that incident plane waves with vanishing en- 
ergy are totally reflected by the barrier. Since the <p{x, k) 
is continuous even at x = ±i?, we consequently see that 



ip(x, ±0) = for x £ 



(27) 



This is just the case without the zero-energy resonance 
0, El • It is worth noting that if we assume that the ini- 
tial wave packet ip(x) is absolutely integrable, Lebesgue's 
dominated convergence theorem leads to 



^(±0) 



<p*(y,±0)ip(y)dy = 0. 



(28) 



Hence the condition for to = 1 can be always sat- 
isfied for this system. This result implies the use of Eq. 



(II 711 with to = 1, which leads to an estimation such that 
ip(x,t) behaves as t~ 3 ^ 2 Q , i.e., 

~ J2 d M^^) ^ (1) H). (29) 

^ ' a— ± 

However, the determination of the asymptotic behavior 
of ip(x,t) may also need a consideration to the concrete 
behavior of the ip( x )- F° r definiteness, we here confine 
ourselves to deriving the behavior like i -5 / 2 , and assume 
that the ip(x) satisfies Eq. (|18fl and has a continuity at 
zero momentum: 



^(«)(+0) =^W(-0) for n = 0,1,2,3. 



TXn) 



(30) 



This allows us to represent both ^W(+0) and ip( n '{— 0) 
by the same symbol ?/»(0). In this case, substitution of 
Eqs. 126|) and (|28() into (|21() leads to a simple expression 
for ^W(±0): 

^\a0) = g*_ (1) (a0) #(0) + 2<5 + ^«(0), (31) 
where g_ (±0) are straightforwardly evaluated as 

5 a) (+0) = 2i J R + 2(cosh2fc bj R)/(iA: fc sinh2fc 6j R), (32) 



<? a) (-0) = -2/(ifc 6 sinh2fc 6 i?). 



(33) 



Equation (|31|l implies the fact that if tp(x) satisfies the 
condition (|22|l for to — 2, i.e., 



V>(0) = and ?A (1) (0) = 0, 



(34) 



then, $W(±0) = and the formula J5§J is no longer 
effective. Note that this does not immediately mean the 
case of to = 2 in the condition (|ll|l realized. Because, 
as is pointed out after Eq. the system's properties 

(|23> with Eqs. (JSDJ and (jS} causes ^ (2) (±0) = 0, even if 
^( 2 )(0) ^ 0. See Eq. Therefore, we should actually 

refer to the case of to = 3 in the condition Ullfl . This 
time, the formula Q17[) (with to = 2) is used again to 
read 

~ \^TI E 9 Mx,<r0) ^ (3) M). (35) 

In this case, ^ 3 )(<t0) = 3 3 l (1) (cr0) ^ (0)+2<5 +i<T ^ 3 ) (0). 

In order to illustrate the above analysis, we consider 
the long-time behavior of the nonescape probability P(t) 



P(t) 



\tp(x,t)\ 2 dx, 



(36) 



instead of ip(x,t) itself. This is the probability of find- 
ing a particle, initially prepared in state ip, in a bounded 
interval / = [a, b] at a later time t. By substituting Eq. 
(|14fl or l|17|) into the definition l|36|) . the asymptotic be- 
havior of P(f) directly reflects that of ip(x,t). We here 
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10° 10 1 10 2 10 3 



t 

FIG. 1: Nonescape probabilities P(t)'s for initial wave pack- 
ets ij} = cj>o, fa, and fa and their asymptotes predicted by Eq. 
1291 or 1351 1 (solid lines). For fa and fa, P{t) shows the well- 
known t~ 3 behavior at long times (long-dashed and short- 
dashed lines), whereas P(t) for fa exhibits another power-law 
behavior like t (a dotted-line). 



restrict ourselves to the three wave packets <j>o(x), fa(x), 
and fa(x), as the initial wave packets ip(x): 

$ m {k) = Nmk ™ e -ao 2 (k-ko) 2 /2-tk XQ ^ (37) 

where m = 0,1,2, ao > 0, fco,£o G K, and N m being 
the normalization constants. Note that the parameters 
ao and xq roughly indicate the width and the location 
of (f) m (x), respectively. These wave packets are rapidly- 
decreasing functions and satisfy the regularity (|30|1 . Both 
fa(x) and fa (x) are chosen to confirm the asymptotic for- 
mula l|29|l . while fa{x) satisfies the assumption (13411 to 
realize the asymptotic formula (|35f) . These wave pack- 
ets do not satisfy the assumption l|18fl . This situation 
however could be taken into account. In fact, if we take 
appropriate parameters ag and xo where the latter satis- 
fies |x | = — x » R, the errors in (/>™^(±0)'s might be 



negligibly small 01 ■ 

Figure ^ shows the time evolution of P(t), and the 
asymptotes predicted by Eq. i|29|) or l|35[) . for the three 
initial wave packets fa(x), fa(x), and faix). In our cal- 
culation, we have chosen a set of parameters ao = 1.0, 
/co = 1.0, and xq = —20.0 for the three initial wave 
packets. Then, every average momentum in these initial 
states is positive. We have also chosen in all these cases 
the potential range R = 1.0, which is much smaller than 
the location |a;o| = 20.0, and height Vq = 16.0, which is 
greater than the expectation value of energy ((f> m , H^fan). 
The interval / — [a, b] for the nonescape probability is set 
around the initial location of the wave packet xq, and we 
set a = —22.0 and b = —18.0. One may recognize three 
different regions in the figure: for small t, all P(t)'s de- 
crease smoothly, and then, they partially revive before 
decreasing again. These regions reflect the motion of a 
wave packet, i.e., it leaves the interval / for the barrier, 
returns to the interval after a collision with the barrier, 
and goes outside through the interval, respectively [l8| . 
It is clearly seen that, in the last region, P(£)'s for ini- 
tial wave packets fa and fa approach to the asymptote 
parallel to the well-known i~ 3 . On the other hand, the 
behavior of P(t) for the initial state fa is in quite agree- 
ment with the asymptote parallel to t~ 5 , other than i~ 3 . 

To summarize, we have considered the finite-range 
potential-systems for one dimension and explicitly char- 
acterized the various power decreasing behaviors of the 
scattered wave-packets at long times, in terms of their po- 
sition and momentum bahavior at an initial time. Our 
results can also cover the free case of Refs. an d that 

of Refs. 0,0 with a slight modification. The power-law 
decrease of the potential systems at long times still has 
not been observed experimentally, however it may exhibit 
an interesting phenomenon, involving a peculiar struc- 
ture where the characteristics of the initial state play a 
crucial role. 

The author would like to thank Professor I. Ohba and 
Professor H. Nakazato for useful and helpful discussions. 
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